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Abstract. We give a rigorous and self-consistent derivation of the elementary braid matrices 
representing the exchanges of adjacent Ising anyons in the two inequivalent representations 
of the Pfaffian quantum Hall states with even and odd number of Majorana fermions. To 
this end we use the distinct operator product expansions of the chiral spin fields in the 
Neveu-Schwarz and Ramond sectors of the two-dimensional Ising conformal field theory. 
We find recursive relations for the generators of the iiTeducible representations of the braid 
group ,'S^2n+l in terms of those for .b^2n, as well as explicit formulas for almost all braid 
matrices for exchanges of Ising anyons. Finally we prove that the braid-group representations 
obtained from the multi-anyon Pfaffian wave functions are completely equivalent to the spinor 
representations of SO(2;j + 2) and give the equivalence matrices explicitly. This actually 
proves that the correlation functions of 2n chiral Ising spin fields a do indeed realize one 
of the two inequivalent spinor representations of the rotation group SO(2n) as conjectured by 
Nayak and Wilczek. 
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1. Introduction 

One fascinating application Uj |2j of the anticipated non-Abelian statistics of the chiral spin 
fields in the critical two-dimensional Ising model has established a remarkable connection 
between the rich-of-exact-results area of the two-dimensional rational conformal field theories 
(CFT) and the new and promising field of topological quantum computation |3J. The localized 
non-AbeUan Ising anyons, which are believed to be reaUzed in the fractional quantum Hall 
state at filling factor v — 5/2, that is most likely described by the Moore-Read (Pfaffian) 
CFT m, turned out to be a useful tool for topologically protected quantum information 
processing lU |2] |5l [51 . Protection against noise and decoherence is obtained by encoding 
quantum information in robust topological characteristics of the strongly correlated electron 
system, such as quasiparticle's fusion channels, while quantum gates are implemented by 
topologically non-trivial operations |[T]|2]|5]|6l, such as braidings of non-Abelian anyons. 

Nayak and Wilczek argued in an insightful paper that the Pfaffian wave functions with 
2n Ising anyons at fixed positions belong to a 2"^ ' dimensional spinor representation of the 
rotation group SO(2«), see Sect. 9 in Ref. Q. However, as explained in Sect. ll.3l below. the 
arguments they presented in support of this claim were incomplete and partly misleading. To 
our knowledge, this claim therefore has been, up to now, only a conjecture. In the present 
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paper we provide a complete proof of this conjecture meeting the requirements of rigor of 
mathematics. This proof contains three steps: 

(i) Derive the elementary generators, b''^'^\ j = 1,2,3, of the braid group ^4, in the 
two inequivalent irreducible representations with positive ("+" in the superscript) and 
negative ("— " in the superscript) fermion parity, directly from the 4-quasihole Pfaffian 
wave functions. 

(ii) Construct recursively the generators 1 < 7 < 2n + 1, of the irreducible 

representations of the braid group ^2n+2 in terms of the generators, B^^"''^\ I < j < 
2n- 1, of .^2«- 

(iii) Find explicitly the equivalence matrices mapping the obtained generators of the 
representations of the braid group ,S§2n+2 to those generating the representations of 
Nayak-Wilczek with the corresponding parity. 

Contrary to the common believe, the first step has not been accomplished by NW in Ref. 17] — 
they have only computed the first line of the 2x2 matrix representing the exchange of the 
anyons with coordinates T71 and T73, which at best could be used to determine the generators 
of the positive-parity representation of ^4, though some ambiguities had to be resolved, see 
Sect. 11.31 below. Later the braid matrix ^^23, as well as the other two generators ^12 and 
Ri4, for the exchanges of 4 Ising anyons have been unambiguously derived in Ref. [jS] by 
careful analytic continuation of the 4-quasiholes Pfaffian wave functions in the representation 
corresponding to even number of fermion fields in the CFT correlation function. 

However, the generators of the negative-parity representation of have never been 
derived before, in the wave-function approach, because the 4-anyon Pfaffian wave functions 
with odd number of Majorana fermions have been unknown. In this paper we give the first, 
to our knowledge, derivation of the generators of the negative-parity representation of ^4 
directly from the Pfaffian wave functions, without computing them explicitly, by using instead 
the short-distance operator product expansions. 

Furthermore, the generators of the two inequivalent representations of the braid group 
^2n+2, spanned by the Pfaffian wave functions realized as CFT correlation functions with 
2« + 2 Ising anyons at fixed positions and even/odd number of Majorana fermions, are not easy 
to obtain. They cannot be derived in analogy with the generators of £§4 because the Pfaffian 
wave functions with 2« + 2 anyons are not known explicitly. Fortunately, it is possible to find 
recursion relations between the generators of ^2n+2 and those of I^2n, by using the fusion 
rules of the Ising anyons. However, when we fuse two Ising anyons, the result could be either 
I or a Majorana fermion \j/ so that in the first case the fusion process maps a representation 
of S§2n+2 with given fermion parity into a representation of ^2« with the same parity, while 
in the second case the fusion process switches to the opposite parity. This subtlety not only 
mixes the representations but also requires that we know both representations of ,'^2n in order 
to construct inductively from them any of the representations of SS2n\-2- 

Meanwhile, the 4-quasihole results of NW have been reproduced in Ref. fS] using the 
universal R matrix in the quantum-group approach for the Ising model. It is worth stressing 
that the arguments of JS) do not prove the NW conjecture because the representation of the 
braid group ^4 in f8^ is defined by the authors in such a way to reproduce the results of NW 
as can be seen from Sects. 5.2 and 5.3 in IS], however, there is no proof that it is the same as 
the braid-group representation obtained by analytic continuation of the multi-anyon Pfaffian 
wave functions as obtained e.g. in f5\. 

In addition, the 4-anyon braid matrices have been convincingly derived in Ref. ||9l for the 
case of the p-wave superconductor, which is known to be related to the Pfaffian state ifTOlfTTl . 
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While the question of the basis orthonormality in the p-wave superconductor has not been 
addressed in ||9l, it has been answered in ifTTl . The analysis of the p-wave superconductor 
is more conclusive about the braid generators, however, the connection to the Pfaffian FQH 
state is more elaborate, because only the large-distance behavior of the weak-pairing phase is 
related to the MR state jTOl . 

The above mentioned confirmations would have been very nice if we had an independent 
proof of the NW conjecture, however, they are still not sufficient to prove this conjecture, 
because they are either defined on purpose to reproduce the NW results or are indirectly 
related to the many-body states of the electron system and the correspondence depends on 
many assumptions. Therefore, it would be useful to have an independent, self-consistent and 
rigorous derivation of the braid matrices directly from the Pfaffian wave functions representing 
the states containing multiple Ising anyons, which are actually used to define the qubits in 

TQc in. 

One more reason for this necessity is that the eventual experiments with the real quantum 
Hall systems would test the properties of the strongly correlated electron state that are encoded 
into the corresponding many-body wave function. Recall that the gauge-invariant quantity in 
the adiabatic transport exchanging Ising quasiholes is the product of the explicit monodromy 
(which can be computed in the CFT or quantum group approach) and the geometrical Berry 
phase lfT2l [T3I [141 which is present only in the wave-function approach. This is a subtle point 
because whether the adiabatic transport of Ising anyons, along complete loops around each 
other, is indeed realized by spinor representations of the rotation group depends on the Berry 
connection of the trial wave functions. It has been argued that the actual holonomy, which is 
the physically observable quantity that we intend to use for topological quantum computation, 
is indeed equal to the monodromy of the Pfaffian wave functions because the Berry connection 
is trivial lfT2l [TTl [T3l [141 . i.e., the only contribution comes from the ubiquitous Gaussian 
factor that is typical for charged particles in magnetic field and this geometrical contribution 
is simply the Aharonov-Bohm phase. Given that the multi-anyon trial wave functions are 
holomorphic, that would certainly be true if they could be proven to be orthonormal |[T3l[T4l . 
While the first attempts [IT2I succeed in generalizing the analogy of (the overlap screening 
of) the Coulomb plasma, at least to the Pfaffian state with two anyons, a recent argument 
about the four-anyon case |[T4l . which could be generalized to more anyons, seems to provide 
convincing evidence that the multi-anyon wave functions obtained in an appropriate CFT basis 
are indeed orthonormal. Therefore, it is now rather plausible that the holonomy of the multi- 
anyon Pfaffian wave functions is precisely given by the monodromy which could be obtained 
by analytic continuation, and this is what we shall use in this paper. Notice also that the 
Landau level mixing, which has important consequences for the physics of the quantum Hall 
state at filling factor v — 5/2, would certainly modify 1 1 3 1 the exchange properties of the Ising 
anyons derived by the monodromies and this effect can only be analyzed in the wave-function 
approach. 

The details in the explicit representation of the braid generators and the differences 
between their distinct realizations become more important when we try to implement various 
quantum gates and to estimate the computational power of the Ising-anyon TQC j5 |. For 
example, it was possible to construct the CNOT gate ID in terms of 7 elementary braidings, 
however this construction could not be generalized for systems with more anyons, i.e., the 
CNOT could not be embedded into systems with more qubits. The precise braid-generators 
analysis is crucial for answering such questions as whether it is possible to implement the 
entire Clifford-gate group purely by braiding or not, see Ref. JS) for the answer 

Outline of the paper: In this paper we give a rigorous and unified derivation of the braid 
matrices, representing the exchanges of Ising anyons, and a proof of the NW conjecture based 
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on the wave-function approach, and this result is independent of the orthonormaHty of the 
CFT blocks used as a basis. Combined with the orthonormality results obtained in Ref. lfT4l 
this implies that the adiabatic transport of Ising anyons could indeed be used for topological 
quantum computation as proposed in I?] H] UT] [141. In Sect. Il.ll we review the subtle issue 
of the chiral fermion parity and its spontaneous breaking in the doubly degenerate Ramond 
sector of the Ising model due to the presence of the Majorana fermion zero mode, following 
Refs. ifTsl [161 and introducing their notation. This is necessary for the formulation of the 
short-distance operator product expansions of the spin fields a± in the Ramond sector, derived 
in ITSI , which is the main tool for the computation of the non-diagonal braid matrices Bj^'*' 
in Sects. [521 and [33] In Sect. 1 1.21 we recall the definitions of the 4-quasihole Pfaffian wave 
functions as correlation functions |[4|[7l in the m(1)x Ising CFT containing four a fields and 
even number of Majorana fermions and introduce the notation for the computational states, as 
well as the encoding of quantum information, following Refs. |[5][6l. In Sect. ll.3l we formulate 
the NW conjecture, which we intend to prove in this paper, and in Sect. 11.41 we announce 
the precise statement about the braid generators that we obtain in the Pfaffian wave-function 
representations of the braid group with positive and negative parity, which is the first step 
towards the proof. Then in Sect. [2] we derive the diagonal braid matrices Bj^'*^ and Bj'*'*' 
by using the short-distance operator product expansion (OPE) in the Neveu-Schwarz (NS) 
sector of the Ising model, representing the elementary 4-anyon exchanges in the positive- 
parity (Sect. 12.1b and negative-parity (Sect. 12.2b wave-function representations. Next, in 
Sect. [3]we derive the non-diagonal braid matrices Bj^'*^, in the wave-function representations 
with positive parity (Sect. 13.2b and negative parity (Sect. 13.3b . by using the short-distance 
OPE, this time in the Ramond (R) sector of the Ising model, which is reviewed in Sect. 13.11 
following Ref. |[T5l . The braid generators B^j''^\ j = 1,2,3, obtained in the representation 
of ,S§4 with positive fermion parity, coincide completely with those derived in Ref. fS) for 
the case when the homotopic condition |J7i2i?34/j?i3i724| < 1 of Ref. fS] is fulfilled. Notice 
that Bj ' obtained here are explicitly different from those of NW. Furthermore, the braid 

generators B^-^' j = 1,2,3, for the representation of 1^4 with negative fermion parity, are 
important new results which have not been obtained before from the 4-quasihole Pfaffian wave 
functions because these wave functions were unknown in the negative-parity representations. 
In the Proposition 1 in Sect. |4| we derive new recursive relations for the projected (2« + 2)- 
anyon exchange generators b^^"^^-'^\ j — 1,...,2« + 1, in terms of those for 2« anyons, 

By^"'*', j — 1, . . . ,2n — 1, and give maximally explicit new formulas for the (projected) braid 
matrices in the Corollary. Finally we prove in the Proposition 2 in Sect.[5]that the braid-group 
representations derived from the multi-anyon Pfaffian wave functions are equivalent to those 
derived in the SO(2« + 2) spinor approach and give the explicit equivalence matrices which 
are also new results. This equivalence makes it completely legitimate to interpret the abstract 
parity in the spinor representations of SO(2« + 2) as the physical fermion parity in the Ising 
model. 

1.1. Double degeneracy of the R-sector and spontaneous breaking of fermion parity 

The Ramond sector (or twisted sector) of the Ising model is defined as the superselection 
sector in which the Majorana fermion field has periodic boundary conditions on the cylinder 
and antiperiodic on the conformal plane 031 [TTl . Because of the periodic boundary conditions 
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of the Majorana fermion in the R-sector, whose Laurent mode expansion in the complex plane 

viz) = ViiZ"'^^^, {Wn. Wm} = 8,,+m.Q 

implies the presence of a fermionic zero mode y/o, the ground state in this sector must 
necessarily be doubly degenerate if the chiral fermion parity is conserved. Indeed, the 
Majorana zero mode i/zij = i//q , (i/o)' = (1/2)1 and the chiral fermion parity operator 7, 
satisfying yVo + VoY~0 and 7^ = I, form a two-dimensional Clifford algebra whose lowest 
dimensional representation is two-dimensional and can be expressed in terms of the Pauli 
matrices lfT6l . Choosing a 7-diagonal basis, the two chiral spin fields of CFT dimension 
1/16 (with positive and negative fermionic parity) intertwining between the vacuum and the 
R-sector's lowest weight state can be written as [131 

|±) = (7±(0)|0), 7(J±7=±(J±. (1) 

The conservation of the fermion parity implies that the two fields a± in Eq. ([T]i must obey 
Abelian fusion rules 

a+xG+ = I, a+x <y_ = Y, c7_ x (t_ = I ^x (j± = (jzp. (2) 

On the other hand, modular invariance requires a single lowest- weight state |[T6|, like in the 
case of the Gliozzi-Scherk-Olive projection in string theory, which is conventionally chosen 
as 

a = ^ — =^ a X a = l+\j/ 

V2 

and consequently obeys non- Abelian statistics. This projection leads to spontaneous breaking 
of the chiral fermion parity, see Ref. lfT6l for a more detailed explanation. Despite the 
seemingly unphysical nature of the chiral spin fields a± with definite fermion parity they 
appear to be very convenient for enumerating different computational states, for labeling the 
fusion paths in the corresponding Bratteli diagrams [3] |6|, as well as for the identification 
of the spinor parity in the representations of SO(2« + 2) with the fermion parity in the Ising 
model. 



1.2. Wave functions for 4 Ising anyons in the positive-parity representation 

The wave function for the Pfaffian fractional quantum Hall state with even number of holes 
(or electrons) at positions zi,... ,Zn and four quasiholes at positions T7i , . . . , J74, can be realized 
as a correlation function, in the m(1) x Ising CFT ||4]|71|5], 

N 

'I'4qh(l?l,I72,I?3,I?4;{z,}) = (rqh(l7l)rqh(l72)rqh(l73)rqh(l74)n'^hole(z;)>(3) 

7=1 

of the field operators corresponding to creation of holes and quasiholes 

V^hoie(z) = vW:e'^^(-^): and v^qh(J7) = (1(77) : e'^T?^^''^, (4) 

respectively, where cy{r]) is the chiral spin field in the Ising model of dimension 1/16 and 
Va(z) is the right-moving Majorana fermion in the chiral Ising model. It can be expressed in 
more explicit form in terms of the Pfaffian wave functions |5J, however, because they will not 
be needed in our fusion-rules approach we will skip them. We shall only use the notation of 
Ref. El m for the two linear independent 4-quasiholes states 

= |0)+ ~ (CJ+(J+CT+(J+), = ~ (a+CT_CT+CT_) (5) 
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and will call them the computational basis in the positive-parity representation (representation 
parity is denoted by the subscript "+"). In the next section we will give more detailed 
expressions for the computational basis states ^ as well as an explanation of the sign 
We can now argue why it is important to derive the braid generators directly from the wave 
functions representing FQH states with 4 Ising anyons. It might indeed be possible to derive 
the braid-group generators from the universal /^-matrix for the quantum group corresponding 
to the Pfaffian FQH state [8 |. However, it can be shown that the Pfaffian wave functions ^ 
with 4 Ising anyons, defined in Ref. 10, are different from the correlation functions of four <y 
fields (which are zero unless eie2 = e$e4 = K and then, see Eq. (6.43) in Ref. ||T5l ) 



1/8 



so that the former functions would have different analytic properties from the latter because of 
the presence of the Majorana fermions. Because the 4-point correlation functions ^ depend 
only on the product of the signs of the fields comprising the first and the second pair we can 
encode information in the topological charge K^e\e2of the first pair and then the topological 
charge of the second pair is fixed to be the same. It is also obvious from Eq. ^ that the order 
of signs in the two pairs is irrelevant, i.e., C7+CJ_ ^ 0+ because e\e2 = eiex = K, so that we 
can always choose the sign of the first a field in each pair to be '+'. 

1.3. The Nayak-Wilczek conjecture 

Nayak and Wilczek conjectured fTl that the elementary matrices representing the exchanges 
of 2n Ising quasiparticles in the Pfaffian fractional quantum Hall state can be interpreted as 
n/2 rotations from SO(2«), i.e., they can be expressed in terms of the gamma matrices 7;'"', 
1 < / < 2«, which satisfy the anticommutation relations of the Clifford algebra 

{^("',7^ }-25,„ l</,;<2„. (7) 

In more detail, the elementary operations for the exchange of the i-th and (/ + l)-th 
quasiparticles could be expressed, in an appropriate basis of 2n-quasiholes Pfaffian wave 
functions ||7j(9i|, as 

D(n) f ^ (") (") A ^'^ fl, (") (n) A .ON 

R) = e'T exp (^--7] V]+'i j ^ ^ (I- 7] VVi ) - (8) 

where 1 < j < 2n — 1 and the second equality follows from the fact that (7;7/+i)^ = —I due 
to the anticommutation relations (|7]i. 

in) 

The 2n matrices 7. have dimension 2" x 2" and can be defined explicitly as follows 

7("' = ai (g) aa (g) • • • (g) ff3 

72"' = (72 (g C73 (g) • • • (g) (73 
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rl"' =I2«-1«)CJ2. (9) 

The "gamma-five" matrix yf\ defined by 

in] ( •\n {n\ {n\ 

7f =(-') rl •■■72„% 

commutes with all matrices ([8]l and therefore 7?^"' cannot change the 7^"' eigenvalues ±1, 
which means that the representation ^ is reducible and the two irreducible components, 
corresponding to eigenvalues ±1, can be obtained by projecting with the two projectors 

Pi"' = 5^,i.e., (Pi:'')'=pr=(pr)'. (10) 

In other words, the generators of the two irreducible spinor representations of the braid group 
SS2n can be obtained by simply projecting (1 < j <1n — X) 

^(«,±) ^ pMj^MpW = ^ (i - 7|"V];\) 4"'. (1 1) 

Eq. ( fTTI ) is what we call the NW conjecture in this paper because the braid generators ( fTTI ) 
have not been derived in Ref. 1|7J from the multi-anyon Pfaffian wave functions. Instead, NW 
say on page 546 in I?]: "We will verify this assertion in the four-quasihole case with our 
explicit wave functions . . . and give an argument in favor of its validity in the 2«-quasihole 
case". To this end they first verify the statement for 4 quasiholes by computing the first row 
of the 2x2 braid matrix representing the exchange /?i3 of the anyons with coordinates rji 
and 773 and then compute R23 from that. However, this kind of derivation of the elementary 
braid matrix /?23 is ambiguous because it is based on the result for the exchange of anyons 1 
and 3 and a braid relation, such as R23 = RuRu^u '^^ ^'23 ~ ^r2'^i3^i2' two results 

are physically different. Formally, this ambiguity appears because the exchange Tji ^ 773 
depends on the homotopy class of the exchange with respect to the second anyon (with 
coordinate 772) and because of the emerging sign ambiguities which have not been fixed 
in a physical way. Next, the generalization argument they mention at the end of Sect. 9 
in Q is as follows: "... imagine bringing 4 quasiholes close together the braiding is 
governed by the OPE and therefore is generated be the transformations we found above in 
the 4-quasihole case". This argument is misleading because NW found the 4-quasihole braid 
generators only for even fermion parity of "the rest", while their generalization argument 
assumes that they could use them also in the negative-parity case which is wrong. The point 
is that when we separate 4 quasiholes the rest 2n — 4 quasiholes could have both positive and 
negative total fermion parity (with equal occurrence in the computational basis). For example, 
consider the 6-anyon computational states given in Eq. (38) in Ref. 15]: following the NW 
"generalization argument", let us concentrate on the first four quasiholes, corresponding to 
coordinates J]i , . . . , 174; the rest of the quasiholes, i.e., the two quasiholes with coordinates 775 
and 775 in this case, have positive total fermion parity in the computational states denoted by 
|00) and |10) in Eq. (38) in f5\, while in the states |01) and |11) it is negative. In the first 
case one could eventually use the four-quasihole results of Ref. Q, while for negative parity 
one needs an inequivalent set of generators which have been missing in Q. Actually the 
braid generators in the negative-parity representations of the braid group £§4 have never been 
known before, because the explicit four-anyon Pfaffian wave functions in the negative-parity 
representation have been unknown. Hence one cannot derive recursively the generators of 
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^2n only from those of the positive-parity representation of M4, so that NW have drawn an 
inference from slight or insufficient evidence. In this paper we are filling this logical gap thus 
turning the insightful NW conjecture into a mathematical theorem. 



1.4. Braid matrices for exchanges of 4 Ising anyons: statement of the result 



One of the important results in this paper is the unified derivation of the braid matrices, B 



(4.±) 



fij^'*^ and B^^'^' , in the two inequivalent representations of the braid group ,^4 corresponding 
to positive and negative fermion parity, directly from the four-anyon Pfaffian wave functions. 
The superscript "4" in the notation for the braid generators expresses that these are generators 
of representations of ^4, while the sign "±" denotes the fermion parity of the corresponding 
representation. In Sections |2] and [3]below we will show that the result for the generators of 
the positive-parity representation of ,^4 



,(4,± 



B 



(4.- 



B 



(4.4 



e'T 
V2 



1 -/ 
-/ 1 



B 



(4r 







(12) 



obtained here by the fusion-rules approach completely coincides with the result of 
Ref. |fT9l obtained by the analytic continuation of the 4-quasiholes wave function when 

|r7l2T734/j7l3'?24| < 1- 

The braid generators B^^'^^ obtained here are explicitly different from those of NW, albeit 

the latter can be obtained by an equivalence transformation generated by (Bg^'^' )^, see Sect. |4] 
below. However, this monodromy transformation makes an observable difference for the 
physical state the topological quantum computer, which has important physical consequences, 
e.g., it implements the Pauli Z gate, and therefore has to be controlled experimentally. 

In addition we shall explicitly derive the generators of the 4-quasihole Pfaffian 
representation of the braid group ^4 with negative fermion parity 



B 



(4,-) 







B 



(4,-) 



e'T 

V5 



1 -i 
-i 1 



B 



(4,-) 



(13) 



which has not been obtained before in the wave-function approach. 

The main idea is to employ the realization of the multi-anyon Pfaffian wave functions as 
CFT correlation functions without using their explicit form. The key point is that the precise 
braid matrices are independent of the distance between the particles being braided because 
they are topological objects. Therefore we could first fuse the particles, which we intend 
to exchange, and then execute braiding by analytic continuation of the relative coordinate. 
For example, the counterclockwise braiding of the quasiparticles with coordinates 771 and T72 
could be executed by the analytic continuation along the circle defined by 



171 + 172 



^! '72 



171+172 



,1712 



0<f < 1. 



2 2 ' '-^ 2 2 

Thus, if we want to exchange the anyons with coordinates Tji and 772 we can first fuse Tji r\2 
inside the CFT correlator, apply the OPE to extract the short-distance singular behavior 
in terms of the relative coordinate 7712 = 77 1 — 772, and then execute braiding by a simple 
permutation 771 ^ 772 plus analytic continuation in 7712, i.e., 

I7i = J72, 172 = J7i, (so that 77(2 = e'''77i2) , 77^ = 77^, for; > 2. 

This leads to crucial simplifications because the (potentially unknown) CFT correlators after 
fusion are independent of 7712 and their explicit form is not needed since the entire non- 
analytic behavior comes from the short-distance prefactors containing 7712. 
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2. Diagonal exchange matrices bJ"^^^ and Bj^^^ fusion rules in the NS sector 

When we exchange the quasiholes with coordinates T71 and 772 or 773 and 774, corresponding 
to the braid generators bJ"*^'*^ and Bj^'*^ respectively, it is obvious from Eq. (O that there are 
even or zero number of <7 fields to the right of the quasiholes being exchanged. Therefore, 
in order to fuse the quasiholes before braiding we can use the OPE of two (7 fields in the 
Neveu-Schwarz sector of the Ising model |T5l[T7l . 

The fusion rules for the a± fields ^ lead to the following short-distance OPE in the NS 
sector ifBlfTTl 

CT±(Z1)CT±(Z2) =^ Zn'^h, CT+(Z1)CT_(Z2) =^ z72^\/^V^fe). (14) 

However there is another contribution to the OPE of the quasihole fields in (|4|i coming from 
the OPE of the Abelian parts of the quasiholes operators 

:e^5fe^(^')::e^5^^<^^>: c zjf : e^^^^^' :, (15) 

— 1/8 

which cancels the factors and from now on we shall skip them (these factors have been 
explicitly shown as T]^^^ in Eq. (9) in Ref. IS). 

It is worth-stressing that the braid matrices Bj^'*^ and B^'*'*^ must necessarily be diagonal 
||5l because the anyons being exchanged are in the NS sector where the chiral fermion parity is 
preserved lfT6l so no coherent superposition of states with different parity is possible. Indeed, 
if we want to exchange 773 with 774 we could first fuse them and then do the braiding. However, 
there are no other a fields to the right of the pair a{T]i)a{T]4) so we have to use the OPE ( fT4b 
in the NS sector and therefore, e.g., the matrix element + (OjBj'*'^' 1 1)+ must be zero. 

In the following subsections we will consider the two cases with positive and negative 
fermion parity separately. 



2.1. Positive-parity representation: 

Using the Abelian CT± fields with definite fermion parity and the fusion-path approach fSjIS to 
label the anyonic states of matter we can write the computational basis Q for 4 Ising anyons 
in the positive-parity representation as follows lfT9l l5] 1211161: 

2N 

|0)+ = ((7+(77i)(J+(772)(T+(773)cJ+(774)nV^(Zi)) 

i=i 

IN 

(cT+(77i)(T_(772)(T+(773)(T-(774)nrfe))- (16) 

i=i 

Recall that quantum information is encoded in the topological charge K of the first pair of 
(7 fields according to the rule |0) ^ (7+ (7+, |1) ^ <7+C7-, while the second pair of a fields 
carries no information - its purpose is to make the total fermion parity in ( fT6b trivial, in order 
for the correlation functions to be nonzero, see Refs. lfT9l l5l 12111611 for more detail. 

(4 ^) 

To compute the braid matrix B\' , representing the exchange of the first two anyons, 
we can first fuse 771 772 and then implement braiding by 7712 — > e"^77i2. The short-distance 
approximation of the two computational basis states are obtained by using the fusion rules 
((11) and O 

2N 

|0)+ ~ (cT+(773)(T+(774)nV^fe)), 
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2^ 

i(V/(772)(T-('73)(T+(r74)nV^(^/))- (17) 

7=1 

Executing the braiding by the analytic continuation 7712 e"^77i2 gives 

2N 2N 
|0)+^ ((T+(773)CT+(rj4)nr(^7))^((^+('73)(T+(T74)nV^fe)) -|0) + , 

1 1)+ ^ /!Iii(^(^2)cT+('73)CT+(7?4) fl V'i^j)) ^ 
Ti ''12 V ^' J 



^(V/(772)c7+(773)(T+(774)nrfe)> ='11>+ (18) 
^ 7=1 

SO that the first braid generator in the positive-parity representation is simply 

"1 0" 
/ 



Similarly, to compute b''^'^^ we first fuse 773 774, using the fusion rules ( fT4b and (fTSl) to 
obtain the short-distance approximation to the computational states, and then braid T734 
s^'rjM to get, completely in the same way (only replacing 7712 with 7734), b''^'^^ = B^^'^\ 

2.2. Negative-parity representation: 

In order to write explicitly the computational basis in the negative-parity representation we 
could introduce one extra Majorana fermion to the right of all a fields, still having even 
number 2N of other Majorana fermions. Thus we define the computational basis for 4 Ising 
anyons in the negative-parity representation 

2N 

|0)_ = (c7+(77i)(7+(772)(T+(773)(J_(774)v/(Z0)n'^(^7)> 

7=1 
2N 

|l)_ = ((J+(77i)(7_(772)(T+(773)CJ+(774)V/(zo)nV^(^/))- (19) 

7=1 

Again quantum information is encoded in the topological charge K of the first pair of a fields 
according to |0) ^ 0'+(7+, |1) ^ C7+C7_, however this time the second pair which fixes the 
total fermion parity of the correlation functions in ( fT9] l, has opposite parity compared to that 
in ( fT6] l, see Refs. 1211 161 for more detail. 

Notice that in general we can insert the extra Majorana fermion between any two 
pairs of anyons. This will define a new basis of computational states in the negative-parity 
representation which is related to ( fT9] l, in which the extra Majorana fermion is to the right 
of all a fields, by a braid transformation that is diagonal, with elements ±1 on the diagonal 
because the Majorana fermion either commutes or anticommutes with any pair of a fields. 

Now we can compute B^' ' by first fusing 771 772 and then taking 7712 e"^J7i2- 
Because the short-distance expansions of c7+(77i)(7±(772) in the NS sector is independent 
of the parity of the other fields in the correlator this braid matrix is the same as for the 

(4 — ) (4 

positive-parity representation, i.e., B\' ~ B\ ' (as matrices because they act on different 
computational bases). 



Ultimate Ising braid-group generators 1 1 

In order to compute Bj^' ' we first fuse 773 774, using the fusion rules (fl4] i and ( fTSl ) to 
obtain the short-distance approximation to the computational states, which gives 

fri — 

|0) c W^((T+(rji)(T+(772)r(rj4)r(zo)nV^fe))' 

'73^14 V ^ j=l 

2N 

|1) ((T+(77,)(T_(772)v/(zo)nV^fe-)>- (20) 

'73 ^^4 

Executing the braiding 7734 — > e"^7734 now gives |0)_ e^^/^jO)- while |1)_ |1)_ so 
that the braid generator for the exchange of the last two a fields in the negative-parity 
representation is 

'/ 0" 



1 



3. Non-diagonal exchange matrices ^2 * ■ Ising fusion rules in the Ramond sector 

Just as in the previous sections we are going to use the fact that the braid matrices for 
coordinate exchanges of two anyons are independent of the distance between them so we 
can simplify the computation by allowing the two anyons to fuse, i.e., letting 7J2 7J3 in this 
case, and reading the exchange phases from the analytic continuation of the singular factors 
containing 7J23. However, when we exchange the quasiholes with coordinates 7J2 and 7J3 there 
is one extra CJ field to the right of the quasiholes being exchanged. Therefore, in order to fuse 
the quasiholes at 7J2 and 7J3 before braiding them we need to use the OPE of two a fields in 
the Ramond sector of the Ising model ITSl [TtI . 

The OPE of two a fields in the Ramond sector of the Ising model is more comphcated 
than Eq. (fl4b because the chiral fermion parity in the R-sector is spontaneously broken lfT6l 
and therefore that OPE might contain more terms. Fortunately this OPE has been explicitly 
derived in Ref. IISJ (see Eq. (6.47) there, in which we identify Ge ='■ (Pe, where e = ± is the 
fermion parity) from the knowledge of the 4-point function computed in Sect. 6 there and 
could be written as follows 

+ (e.e2)^w{V^) iSe,.e2 \-e)+ | + • • (21) 

Recall that in the notation of Ref. ifTsl the ket-vector \e) is defined as the lowest- weight state 
in the R sector with fermion parity e, i.e., \e) :~ (7(,(0)|0). We shall use Eq. (l2Tl i in the 
next Subsections to derive the short-distance approximation of the computational state's wave 
functions in the bases {|0) + ,|1) + } and {|0)_,|1)_}. 

3.1. Short-distance OPE of the computational basis wave functions 

In order to simplify the analysis of the fusion process 7J2 7?3 we shall denote the R-sector 
states entering the 4-qh wave functions ( fTSI l as 

2A' 

l±):-c^±('?4)nV^(^^')|0), 
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and recall that by construction the number 2N of Majorana fermions is even. Let us now apply 
the R-sector's OPE i2T[ . for e = +, to obtain the short-distance expansion of the 4-anyon 
computational basis vector |0)+ (in the positive-parity representation) defined in Eq. ( fT6] l 

CT+(T72)cT+(T73)|+> ^ {'+> + ( + • + ) V^V/(V^7^)h 

Then, multiplying from the left by (0| CJ+ (tji ) we obtain 

|0>+„=^„ ^^((0|ct+(77i)|+> + a/1^(0|(T+(770|->|. (22) 



''2^13 Vlrili^ I V 2 

— 1/8 

Notice that the overall phase factor TJ23 in the above OPE is exactly canceled by the 
additional inverse factor coming from the OPE of the Abelian part of the Ising anyons, i.e., 
from Eq. ( fTSl ) for zi = J72 and zi ~ fjs, and we shall remove it from all expressions below. 
Thus, recovering the detailed notation, we get the OPE of the first computational-basis state 
to be 

I ( 2N 



|0) + ^^^^^ ^ <{ (0| a+ (771 )(T+ (774) n V^fe) |0) 



^{0\a+{mMVmT3)(y-iri4)Uw{zjm \ . (23) 



;=i 

Similarly, for the computational-basis state we obtain the short-distance expansion 

1 f 2^ 
|1)+ C — (0|(T+(T7i)(T+(774)nrfe)|0) 

m^Tli V2 [ y=l 

- y^(0|cT+(T7,)WV^)c^-('?4)flr(z,)|0>|. (24) 
Adding and subtracting Eqs. ( |23] | and ( |24] | we obtain 

M±^^ ^^^^^ -Ly^(0|(T+(770r(V^)^-('74) n rfe 

Eq. ( l25T l will be our starting point for the derivation of the braid matrices in the next subsection 
because it expresses the correlation functions on the RHS in terms of the computational basis 
in the LHS in the short-distance limit. 

3.2. Braiding r\2 with TJ3 in the positive-parity representation 

(4 +) 

The braiding transformation is represented by the coordinate exchange 

J72->r73, ni^ni. so that 7723 ^ e"^ 7723- 
Applying the coordinate exchange over |0)+ and making analytic continuation in T723 we get 

if 2^ 

^2 io)%,s,,7f voi^+(^i)^+(^4)nrfe)io) 



e''^T723 



2N 



|(^+(r7i)V^(%/^)(^-(J74)nrfe)|0)| . (26) 
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Now, using Ve^ = /, we can substitute the correlation functions appearing in the RHS of 
Eq. ( |26] | with the expressions in the LHS of Eq. dZSl l to get 



|0>- 



II) 



,.|0)+-|1> 



m^m 2 2 
Repeating the same procedure for the computational-basis state 1 1 ) + we obtain from Eq. (|24] | 



B 



(4,4 



ID 



|0)+ + |1)+ .|0)+-|1)_ 



SO that the braid matrix in the basis {|0) + , |1)+} is 



B 



(4.- 



1 


'1 + 


• 1 - /■ 


e'f 


■ 1 


2 


1 - 


■ 1 + / 


' V2 


-/ 1 



(27) 



The braid matrix obtained in Eq. ( |27l ) completely coincides with the one obtained in Ref. Q 
for |i7i2T?34/''7i3''724| < 1 ; where it was denoted as R^23- 

3.3. Braiding r\2 with 773 in the negative-parity representation 

Again the computational basis in the negative -parity representation is given by Eq. ( fT9] l where 
2N is even. Let us now denote 

2N 

(j±(T74)v/(zo)nV^fe)|0> =:|T> 
i=i 

and apply the OPE ( 1211 1 to obtain the short-distance version of |0)_, i.e., 



o-+('?2)o-+(r?3)|- 



1 

V2 



w{\/n7Jh) 



and similarly, for the computational-basis state |1)_ we can use the OPE (I2TI 1 in the form 



c^-('?2)o-+(r?3)|-> 



V2 



1+) 



w{\/n7Jh)\-)\ , 



to obtain (after adding and subtracting the results for the two short-distance approximations) 



|0)- 



II) 



-y 2N 

-= (0| (J+ (Tji (TJ4) v/(zo) n V^fo) |0) 

V2 j=i 



\0)- 



m^m V2 V 2 



2N 



Because the above equation has the same 7723-short-distance structure like Eq. ([25), executing 
tu^ u,„;^;„„ , ^i7t„ ^„t,;^ d(4> ) _ Z?!'*'^' (as matrices). Thus we 



the braiding 7723 s"^Ti23 produces the same matrix 



conclude that B^' ' is indeed given by Eq. ( fT3] ). 



4. Braid generators for exchanges of 2n + 2 Ising anyons: wave-function approacli 

To summarize our results for the exchanges of 4 Ising anyons, we note that the generators of 
the negative-parity representation ( flJl l of ^4 completely coincide with those obtained directly 
from the / matrices for SO(4) in Sect. Ill of Ref. |I2T1, i.e.. 



B 



(4-) ^^(2,-) 



= 1,2,3. 



Ultimate Ising braid-group generators 



14 



In contrast, the second generator of the positive-parity representation of IJ§4 is different from 
that obtained in the 7 matrix approach. Nevertheless, it is easy to see that the positive-parity 
representation (fTSI i. obtained from the 4-quasiholes Pfaffian wave functions, is completely 

(2 +) 

equivalent to the positive-parity representation in the /-matrix approach Rj ' (see Eqs. (9) 
and(10)inRef. ||2D), i e., 

Bf+)=Z/?f+'z, ; = 1,2,3, Z=(<±))' = 

and the equivalence transformation is explicitly given by the Pauli matrix Z. This is so because 

(4 +) (4 +) 

the Z matrix commutes with the diagonal B\ ' and By , however, changes the signs of the 

off-diagonal elements of ^2^'^'. It can also be directly seen that the two representations ( fT2] i 
and (flJl i of ^4 are inequivalent, see e.g.. Sect. Ill in Ref. 1211 . Thus we conclude that the 
representations ( fT2] i and (fTsT i of ^4 are equivalent to the spinor representations of SO(4) with 
the corresponding parity llTllSTI. 

In this section, we shall generalize this result to the braid representations for 2n + 2 Ising 
anyons. Our strategy to compute the braid matrices b^J"^'^''^\ describing the exchanges of 
2« + 2 anyons in the Ising representation of the braid group ,'^2n+2^ would be to fuse some 
pair of (7 fields, representing one of the qubits in our n qubit system, which has the effect 
of projecting out this qubit. The resulting states after fusion will belong to one of the two 
representations of with positive or negative parity so that we can express the braid ma- 
trices B^" ' recursively in terms of B^ . More precisely, we shall prove the following 
recurrence relations: 

Proposition 1: the 2" x 2" dimensional matrices (1 < y < 2«+ 1) representing 

the generators of the braid group ^211+2 in the computational bases ( |36] | can be expressed 
recursively in terms of the braid matrices B^^"'^\ (1 < y < 2n— 1) generating the Ising 
representation of ^2n as follows: 

(i) ^(2„+2,+)^^(2«+2,-) l<;<2n (28) 

(ii) 5(2"+2.±) ^ g(2n.±) ^ ^^j. 1 < j < 2„ - 3 (29) 

(iii) Bf"+''±'=Bf_"2±'©Bff' for 3<7<2«+l (30) 

Proof: We shall prove this proposition by induction with a base n = 2. To this end we shall 
first explicitly prove statements (i)-(iii) for n = 2. The braid generators B^^'*^ are 4x4 
dimensional matrices defined in the computational basis for 6 anyons (corresponding to 2 
qubits, encoded in the first two pairs of a fields, plus one extra inert pair, formed by the last 
two a fields). For positive parity this basis can be written as |l6l 

|00)+ = {a+o+a+a+a+o+}, |01)+ = {a+o+o+O-a+a-), 

|10)+ = (C7+C7_C7+C7+C7+C7_), 1 11)+ = ((7+<7_Cr+Cr_C7+C7+), (31) 

where we skipped for simpUcity the product of the even number of Majorana fermions as well 
as the coordinates rjj of the fields (J±{rij). In order to find the braid matrices B^'^^ or B^,^'^', 
exchanging rji ^ 772 or 772 TJ3, respectively, we can first fuse 775 — > 776. The results after 
fusion are computational states from the positive- or negative-parity representations of ^4, 
i.e., 

|00)+ ^ ((J+(J+(7+(J+) = |0) + , |01)+ ^ ((j+(7+(j+(j_v/) = |0)- 

|10)+ ^ (C7+(J_C7+(J+V/) = |1)-, |11)+ ^ ((J+CJ_(J+CJ_) = |1) + . 
'75^176 175^176 



1 
-1 
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Let us first compute Using the above fusion results for the computational basis, as well 

as Eqs. ( fT2] i and ( fTSl l for B^^'^\ we obtain 



Bf•+'|00)^ 
Bf'+)|01)H 
Bf'+)|10)H 
<+)|11)h 



(4,4 



SO that B 



= diag(l,l,;,/) =B^ 

44,.) 



:Bf-)|0) 
:fir-)|l) 



|0)+ = |0) + 
|0)- 
'11)- 
'11)- 



|00) + 

|oi)+ 

i/|10)H 

./|11)h 



(8)l2. Next, in the same way we compute B^f'^' by 



using Eqs. (HDi and ([T3ll forB^'* '^^ i.e., we have B^°'^^|00)_ 



BY'^'\0)+ so that 



|00)h 



e'T 



(|0)+-'11)- 



e'T 

71 



(|00)+-/|11)+). 



Similarly, we have Bj 1^1)- 



e'T 



B^"*' so that 



e'T 



|01)+ _(|0)_ -/|1)_) c — (|01)+ -/|10)+) 



V2 



Continuing in this way with the states 1 10)+ and 1 1 1) 

1 

,(6,+) _ ' 
^/2 



V2' 
we find 



B 







1 

-/ 
l-i 



(32) 



Next, in order to compute the rest of the braid generators we can instead fuse the first two 
Ising anyons. This projects out the first qubit so that the computational basis becomes 



|00)h 
|10)h 



(c7+a+(j+(j+) = |0)+, |01)_ 

(V/(T+(T+CJ+(T-) = |0)_, 111)- 



{a+a-a+a-) = |1)- 



(V/C7+CT-(T+(T+) = |1)-. 



(6r 



Consider, e.g., the braid matrix Bi' . It is obvious that 



B 



(6,+) 
3 

3 



|00)h 



Br^^ 1 10)- 



iB 



(4.- 



10)- 



Bf^ 
Bf- 



'101)- 
'|11)^ 



B 



(4,- 



B 



(4,- 



'|1) + , 

'ID-, 



so that B 



(6,4 



= B 



(4.4 



)B 



(4,- 



I2 ® B 



(4,4 



because B 



B 



(4,- 



(33) 

Here we used 



the sign © to denote the direct sum of matrices. Completely in the same way we find 



B 



(6,- 



B! 



(4,4 



)B 



(4,- 



h<»B^2'' 



because B 



(4,- 



= B 



(4.- 



andB 



(6.4 



B 



(4,- 



)B 



(4,- 



diag(l,/,/, 1). Notice that the last generator Bj^'^' is not a tensor product of I2 and B3 



?(4r 



because Bg"*' ' 7^ Bj 

Next we have to repeat the above computation of the braid generators By 
negative-parity representation. The computational basis is now given by [61 



,(4,-H 



(6,- 



|00)_ = (a+a+a+a+a+o-) , |01)_ 

|10)- = (cT+(7_C7+(7+C7+(7+), |11)_ 



(c7+<T+cr+cr_CJ+CJ+), 

( (7+ CT_ C7+ (7_ (7+ (7- ) , 



in the 



(34) 
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because the total parity of the a fields must be negative (and this is compensated by the odd 
number of Majorana fermions inside the CFT correlators which are omitted again). It is not 
difficult to see that the results for B^^' ' are very similar to those for Z?^^'^' just in each step 

all + are replaced by — and vice versa. For example B^"' ' = diag(l, 1,/,/) = ' (E)l2, 
^(6,-) ^ ^(6,+) ^ consequence of 4^-' = = ^i''"' ©^f'^' = h^^B^^''^ , 

etc. The only difference is in the last generator where Bj^' ' = B^'' ' © Bg^'^' = diag(/ , 1,1,/) 
cannot be written as a tensor product. 

The results for the braid matrices B^.^'*'' can be summarized as follows: B^^'*^ are equal 
and given by ( |32] | and the others are given explicitly by 

Bf-±) =Bf'±)®l2, Bp) =Bf'±)®B(/'^' =l2®Bf'±\ 
Bf±)=l2®B(^'±\ Bf±)=B(^'±)eB(^'^\ (35) 

where B^-^'*^ are defined in Eqs. (fT2] | and il3[ . This proves statements (ii) and (iii) for the 
case « — 2, which is our induction base. Now it is easy to see that in addition these braid 
matrices satisfy 

^(6,+)_^(6,-) 1<;<4, 

J J — J — 7 

which proves the statement (i) for the base n = 2. 

It can also be seen that the positive-parity representation of the braid group obtained 
here is completely equivalent to the one derived earlier in Refs. |fT9l l5l l22l and the equivalence 
is estabhshed by the braid matrix U = b'"^"^'^ Bf'^'^ B'f'^'^ Bf'^"^ representing the exchange of 
the pairs (773 , 774) and (175 , 776) (recall that in the representation of Refs. |fT9l |5] |22l the inert 
pair was G{rji)G{rj4) while here the inert pair is a{rjs)G{ri(,))- 

Induction step: Let us assume that the statements (i)-(iii) are fulfilled for the matrices Bf"^\ 
We must first specify the basis of computational states for 2n + 2 anyons in which the braid 
matrices are represented. The general scheme for representing n qubits in terms of 2n + 2 
Ising anyons could be described as follows 11211 16] . We group the 2n + 2 fields a into n+l 
pairs and encode information into the first n pairs: the state of the i-th qubit is |0) if the i-th 
pair of (7 fields is cj+(t72,-i)c7+(772i ) (i-C-, it fuses to the channel of I) or |1) if the i-th pair 
is 0'+('72(-i)o'-(J72i) (i-C-, it fuses to the channel of \ir). The last pair C7+(T72n+i)(7t.(T72„+2) 
contains no information because its state c is determined by the requirement to have a non- 
zero CFT correlator, i.e., c = H^^i c,. Thus, the computational states in the positive/negative- 
parity representation of our «-qubit system are defined as CFT correlation functions of the 
(2n + 2) non-Abelian a fields and an even/odd number of Majorana fermions 

N 

|ci,...,c„)± = {a+Oc, ■■■a+ac„(y+acY[w{zj))- (36) 

j=i 

The parity of the representation is denoted by the subscript of the computational basis states: 
it is '+' for positive parity (corresponding to even number of Majorana fermions) and ' — ' 
for negative parity (corresponding to odd number of Majorana fermions). In other words, 
Cj = + corresponds to the state |0) of the j-th qubit, while cj = — corresponds to the state 1 1 ) . 
Following our strategy we can first fuse the fields G{rj2„-\) and (7(772,,) corresponding to the 
last qubit which has the effect of projecting out the last qubit, i.e., 

I , J |ci,C2,...,c„_i)± ifc„ = + 

mn-l^nin { ±|Cl,C2,...,C„_l)^ itc„ = - 
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(The sign ± in front of the (« — l)-qubit computational state for c„ = — coincides with the 
eigenvalue of the braid generator ^j^"^"}"', however, it is unimportant for our purposes and 
we will skip it below.) This means that the computational states after projection will be 
organized in pairs, such as |0, 1,0. . . ,0, 1, 1)+, |0, 1,0. . . ,0, 1, having exactly the same 

^ V ' V ' 

n-l H-1 

State of the (« — 1 )-qubit system, however with opposite parity, on which the (2« + 2)-particle 
exchange would act by B^.^"'^' and B^^"' \ Given that we do not touch the last 4 

anyons, corresponding to the last qubit and the inert pair, and provided that the braid matrices 
acting trivially on the last 4 anyons in the two representations are the same B^^"'^' = B^^"' \ 
(1 < 7 < 2« — 2) because of the inductive step (i), we arrive at Eq. ( |29] |. 

On the other hand we can fuse instead the first two a fields corresponding to projecting 
out the first qubit. Then, in the first half of the computational states, containing <T+(7+ as a 
first pair, the result of fusion is I so that the remaining CFT correlation function describes 
a (n — l)-qubit computational state with 2n anyons and the same parity. In the second half 
of the computational states, containing a+G as a first pair, the result of fusion is y/ so that 
the remaining CFT correlation function describes a (« — 1 )-qubit computational state with 2n 
anyons however with the opposite parity compared to the original one, i.e., 

T],^J72 [ -|c2,...,c„_i,c„)^ ltci=- 

The minus sign multiplying the (n — 1 ) -qubit state when ci = — 1 , which is totally unimportant 
here because the braid generators which we want to compute act linearly, comes from the fact 
that after fusing the first two anyons in this case we get one Majorana fermion on the left of 
all remaining a fields, which we have to move, according to our convention, all the way to the 
right of them as it is in the definition ( l36b of the computational states with negative parity. This 
produces one minus sign for each pair of a's which is in the state 1 1) but the total sign for this 
move is always ' — ' . Next, executing exchanges on the remaining 2n anyons, that do not touch 
the first qubit, we immediately find the recurrence relations (|30] | (notice the shift j ^ y — 2 in 
the indices of the braid matrices due to renaming of the remaining anyons coordinates TjJ = 
Tlj-2 for 3 < j <2n + 2). The above mentioned extra minus sign does not change anything 

because the braidings act linearly. To illustrate this consider, e.g., the action of Z?g on the 
states |000)+ and 1 100)+, which after projecting the first qubit will go to±|00)±. We have for 

the first state fif +'|000)+ ~ fif (|00)+) = 2^ (|00)+ - /|01)+) = £^ (|000)+ ^ /|001>+), 
while for the second one 



C"^|ioo)+ ^ Bf-'(-|oo)_)^-:-=(|oo)_-/|oi)_) ^ 

-^(-|100)+ + /|101)+) = ^(|100)+-/|101)+). 

Finally we must prove statement ( |28] | for the matrices Indeed, we have 

assumed that B^^"' ' = f^j. i < j <2n — 2, which is the inductive step (i) for bJ'"''^\ 

Then we can first consider the case when 1 < y < 2n — 3 and use ( |29] l to find 

^(2„+2.-) ^ j^n.-) ^ ^ ^(2„.+) ^ ^ j^[2n+2.,+)^ 1 < j < 2« - 3. 

For the rest of the braid matrices we can use (iii) to prove that for 3 < j <2n 

g{2n+2.-) ^ ^(2n,-) ^^(2n,+) ^ ^(2n,+) ^ ^(2«,-) ^ ^(2n+2,+) 



e'7 
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for 1 < / < 2n — 2 where / — j — 2. This completes the proof of 



because B*-.^"'^' 
the Proposition. 



Corollary: The recurrence relations (i)-(iii) in the Proposition allow for the following exphcit 

It be useful: 

12, for 1 < < n, (37) 



representation for most of the braid matrices which might be useful: 



Bg."_+2.±) ^ J2,^...^]l2( 



j-1 



1 
i 



■' ^ , ' V2 

j-i 



for n>2 and 1 < y < n — 1 , as well as 



J2n+2.±) 
'2n 



n-l 



e'T 

7! 







n-J 




- 1 













1 


— ; 








— / 


1 





.— i 








1 . 


" 1 


— / 






— ; 


1 


7 





(g)l2(8)---(8)l2, (38) 



(39) 



plus only one more (non-trivial) recursive relation for the last (diagonal) generator 



B 



(2n+2,±) 
2n+l 



B 



(2«,±)^^(2«,=f) 



2n-l 



2n-l 



(40) 



with a base b\ ' , given in Eqs. (fT2] | and ( fT3] l. Notice that braid matrices ( l40l i are different 
in the two representations with even or odd number of Majorana fermions and cannot be 
expressed as tensor products of braid operators form the braid groups, such as .^^n or ^4, for 
smaller number of anyons. To illustrate the derivation of the explicit formulas ( |37] |. ( |38] ) and 



( [39] ) consider for example the matrix B^ 
have from Eq 



(2n+2.±) 



(i.e., y = 2 in Eq. (l38Tl): on the one hand we 



B 



(2n+2,±) 



(2«.±) 



) I,,,-, 



2"-3 , 



where we have used that we can add tensor factors of I2 to the right, reducing at the same time 
the value of 2«, until 4 < 2« — 3, i.e., until 2n > 8, which gives rise to (2n — 8)/2 + 1 factors 
and on the other hand, using Eq. ( l30l l. we have 



j(2«+2.±) 



ii2«- 



}B 



(6,±) 



(g)l2«- 



because B 



(8,±) 



,(6,±) 



:I2® 4^'"^' (note that 4^'^ 



Eqs. ( I30I 1 for y = 2n and ( l28l l we can verify Eq. ( |39] | 



bI^' '). Similarly, combining 



B 



(2;j+2.±) 



B 



(2«,±)^^(2«+2,t) 



I.. 



)B 



(4,±) 



-■2,1 --"2n-2 ^"2n-2 —'^2V.'^2n-2 —'^2"-i-^"2 

Equations (IJTT i. (l38T l. (|39] l and (l40l i give the most exphcit expressions for the generators 

^(2«+2,±) ^^jj^g ^^^.^ 

group •'^2n+2 in the two Ising-model representations with opposite par- 
ity. These equations also allow us to express the braid matrices, representing the exchanges of 
Ising anyons in the multi-anyon Pfaffian wave functions, in terms of the universal R matrix for 
the Ising model |l8]|22l (or, equivalently, for the su{2)2 Wess-Zumino-Witten model). Note 
the crucial role of the projectors to states with definite parity leading to topological entangle- 
ment 122], i.e., to the fact that not all braid generators are expressible as tensor products of 
braid matrices with smaller dimensions and the unit matrix I2. 
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5. Proof of the Nayak-Wilczek conjecture 

Proposition 2: The representations of the braid group ■S§2n+2 with positive or negative 
fermion parity, constructed from the multi-anyon Pfaffian wave functions, are equivalent to the 
spinor representations, with the corresponding parity, of ^2n+2 constructed from the spinor 
representations of SO(2« + 2). In more detail, the generators 5^.^"+^'*' of ^2n+2 in the wave- 
function representations can be expressed in terms of the generators 7?*"^' *^ in the SO(2«-|-2) 
representations as follows 

^(2„+2,±) ^ |^^(2„+2,±)^ -\(n+l,±)^(2n+2,±)^ 1 < J < 2« + 1 , (41) 

where the equivalence matrices for the positive/negative parity (+/— ) are given explicitly by 
the product of all diagonal generators 

C(2„+2,±) ^ Y\ R^]+\-^^ = n B2j'-i^'^^ ■ (42) 

Proof: First of all it is easy to see that the diagonal generators in the wave-function and spinor 
representations of ^2«+2 coincide, i.e., 

„(2«+2,±) _ p(n+l.±) , ^ ; <- „ I 1 

Indeed, as we can see from 1211 . the diagonal matrices fig"^^'*' with indices 1 < 2y — 1 < 
2n — 1, given explicitly in Eq. dSTl l above, are completely identical to the diagonal matrices 
^g-i'*^ for 1 < 27 - 1 < 2n - 1, given explicitly in Eq. (26) in Ref. ||2T|. In addition, the last 
diagonal matrices are equal because they satisfy the same recurrence relations 

„(2«+2,±) _ „(2n,±) „(2«,=F) _ p("+l-±) _ n(n.±) rr, d("'=F) 

with exactly the same bases, b'-j'^^ = R^-'^l Bee ause the equivalence matrices ( |42] | are 

diagonal by construction, the matrices fig"^^'*' and ^g^J'*' trivially satisfy Eq. (HTI) . and 
therefore we only need to consider the non-diagonal matrices. The non-diagonal matrices 
in the SO(2n + 2) representation r["+^'^^ = ^^j^ expressed El] as 

projections of the unprojected matrices (c7i and (72 below denote the Pauli matrices) 

/?g.+ '' = Il2(8) — (8)l2 g)--^ (I4 - i02 «) C72) <g) l2'g)- - -(8)l2 . 

On the other hand, Eqs. ( [38] l and ( |39] l suggest that the braid generators fig"^^'*' can be 
expressed in a similar way as projections Bg"^^'*' = p^"+^)gg."+^)p£'+') (with the same 
projectors as forTjg^''^^ given in Ref. 1211 ). of the unprojected matrices 



=l2«)---Ol20^(l4-iai«)ai)«)l2«)---®l2, (43) 



i.e., if we define the unprojected matrices fig"^^' as in Eq. (|43l l and apply the projectors f^'^'' 

as described in flT] then the projected matrices will completely coincide with Zjg"^^'*' as 
given in Eqs. (|38T l and (39[ . Next, we can directly prove that the unprojected matrices are 
related by 

^(2„+2) ^ J^^(2„+2)^ ^(«+i)c(2„+2)^ 2 < 2; < 2«, (44) 
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where the unprojected conjugation matrix is 

This is simply because of the identity 

Now, projecting both sides of (l44l i with the projectors taking into account that 

C{2n+2)p^'+i) ^ p(n+i)(.(2n+l) ' ^ p^^'+i', we obtain Eq. (ED where the 

projected equivalence matrix is equal to the product of all diagonal projected braid matrices 
and coincides with Eq. ( |42] |. which completes the proof of Proposition 2 . 

Propositions 1 and 2 ultimately prove that the Pfaffian correlation functions with 2n + 2 
non-Abelian quasiholes at fixed positions indeed belong to one of the two inequivalent 
representations of the braid group ^2n+i whose generators b'^j'"^^'^'' can be expressed as 
n/2 rotations in terms of the SO(2« + 2) y-matrices and identifies the parity in the spinor 
representations with the fermion parity in the Ising model 1711211161. 

6. Conclusions 

In this paper we have consistently derived the braid matrices representing the exchanges of 4 
non-Abelian Ising anyons in both representations with positive and negative fermion parity. To 
this end we have used the fact that the braid matrices are independent of the distance between 
the braided particles, as well as the fusion rules for the Ising anyons in both Neveu-Schwarz 
and Ramond superselection sectors of the Ising model. In addition we found recurrence 
relations for the braid matrices b^J'"^'^''^^ for the exchanges of 2« + 2 Ising anyons as well 
as explicit formulas for most of the braid generators in the representations with both parity. 
Finally, we have proven that the braid matrices derived from the multi-anyon Pfaffian wave 
functions are completely equivalent to the braid generators derived in the SO(2« + 2) spinor 
approach I?] |2T] and have given explicitly the matrices establishing the equivalence in both 
representations. 
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